perm filename DEG6.IN[225,JMC] blob sn#005373 filedate 1971-03-05 generic text, type T, neo UTF8
00100	MATRIX MZ(5,5),M12(5,5),M34(5,5),M23(5,5),M45(5,5),M56(5,5);
00200	
00300	FOR I←1 STEP 1 UNTIL 5 DO
00400	FOR J←1 STEP 1 UNTIL 5 DO
00500	MZ(I,J)←0;
00600	
00700	M12←M34←M23←M45←M56←MZ;
00800	
00900	M56(1,1)←C1;M56(2,2)←C2;M56(3,3)←C3;M56(4,4)←C4;M56(5,5)←C5;
01000	
01100	M45←MAT((X11,X12,X13,X14,X15),(X21,X22,X23,X24,X25),
01200	(X31,X32,X33,X34,X35),(X41,X42,X43,X44,X45),(X51,X52,X53,X54,X55));
01300	
01400	M12(4,4)←M12(5,5)←M23(4,4)←M23(5,5)←M34(4,4)←M34(5,5)←1;
01500	
01600	FOR I←1:3 DO FOR J←1:3 DO BEGIN
01700	M12(I,J)←B12(I,J);M23(I,J)←B23(I,J);M34(I,J)←B34(I,J) END;
01800	
01900	OFF ECHO;
02000	
02100	END;
02200	
02300	C1←C2←C3←C4←1; C5←-1;
02400	
02500	M45←MAT((3/4,-1/2,0,-3/4,1),(-1/4,1/2,0,-3/4,1),
02600	(0,0,1,0,0),(-1/4,-1/2,0,1/4,1),(1/8,1/4,0,3/8,1/2));